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Abstract. We can establish a new picture, the perfect fluid sQGP core and 
the dissipative hadronic corona, of the space-time evolution of produced matter 
in relativistic heavy ion collisions at RHIC. It is also shown that the picture 
works well also in the forward rapidity region through an analysis based on a 
new class of the hydro-kinetic model and that this is a manifestation of rapid 
increase of entropy density in the vicinity of QCD critical temperature, namely 
deconfinement . 
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1. Introduction 

Recently, physicists at Brookhaven National Laboratory made an announcement in 
the American Physical Society annual meeting at Florida, USA, April 18, 2005, that 
"RHIC serves the perfect liquid" [|Tj. The agreement of hydrodynamic predictions 
[E] of integrated and differential elliptic flow and radial flow patterns with Au-|-Au 
data at RHIC energies [El E El EI is one of the main lines of the announcement. 
We first study the sensitivity of this conclusion to different hydrodynamic assump- 
tions in the hadron phase. It is found that an assumption of chemical equilibrium 
with neglecting viscosity in the hadron phase in hydrodynamic simulations causes 
accidental reproduction of transverse momentum spectra and differential elliptic 
flow data in a way that chemical equilibrium imitates a sort of dissipation. From a 
systematic comparison of hydrodynamic results with the experimental data, dissi- 
pative effects are found to be mandatory in the hadron phase. Therefore, what is 
discovered at RHIC is not only the perfect fluidity of the strongly coupled quark 
gluon plasma (sQGP) core but also its dissipative hadronic corona surrounding the 
core of the sQGP. Along the lines of these studies, we develop a hybrid dynamical 
model in which a fully three-dimensional hydrodynamic description of the QGP 
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phase is followed by a kinetic description of the hadron phase [ . We will show 
rapidity dependence of elliptic flow from this hybrid model supports the above pic- 
ture in Sec. 3. Finally, we will argue in Sec. 4 that this picture is a manifestation of 
deconfinement transition, namely, a rapid increase of entropy density in the vicinity 
of the QCD critical temperature as lattice QCD simulations have been predicted [ 

2. SQGP CORE AND DISSIPATIVE HADRONIC CORONA 

A perfect fluid in the QGP phase is assumed in most hydrodynamic simulations. 
While one can find various assumptions in the hadron phase [|2|, e.g. (1) an ideal 
fluid in chemical equilibrium (CE), (2) an ideal fiuid in which particle ratios are 
fixed (or in partial chemical equilibrium, PCE), or (3) a non-equilibrium resonance 
gas via hadronic cascade models (HC). Hydrodynamic results are compared with 
the current differential elliptic flow data, V2{pt), in Fig. 20 in Ref. with putting 
an emphasis on the difference of these assumptions in the hadron phase. The classes 
CE and HC reproduce the data for pions and protons well. This is considered to 
be one of the evidence of success of an ideal QGP fluid at RHIC. Contrary to its 
success, results from the second class, PCE, deviate from the other hydrodynamic 
results and experimental data though the class PCE as an ideal fluid describes 
chemical compositions of the hadronic matter in a more realistic way than the class 
CE does. In order to claim the discovery of perfect fluidity from the agreement 
of hydrodynamic results with V2{pt) data, we need to understand the difference 
among hydrodynamic results and the deviation from data within the PCE assump- 
tion. The difference of assumptions in the hadron phases reduces to the difference 
of the final slope of differential elliptic ffow dv2{pT)/dpT- ^2 is roughly propor- 
tional to pt in low pt region for pions. In such a case, the slope of V2(pt) can be 
approximated by V2/{pt) since one can easily show dv2{pT) / dpr — V2/{pt) when 
V2{pt) is exactly proportional to px- Integrated V2 is generated in the early stage of 
collisions. Whereas differential V2 can be sensitive to the late hadronic stage since 
dv2{pT) / dpT ~ V2I {pt) indicates interplay between elliptic flow (integrated V2) and 
radial flow (mean transverse momentum (pt))- 

In Fig. ^ thermal freezeout temperature T*"^ dependence of (pt) for pions 
including contribution from resonance decays are shown from hydrodynamic simu- 
lations [HII. Here the impact parameter 6 = 5 fm and collision energy ^snn = 200 
GeV are assumed in the hydrodynamic simulations, (px) for pions in the chemi- 
cally frozen hadronic fluid decreases with decreasing T'^ (increasing proper time t). 
This is due to longitudinal pdV work done by fluid elements. Whereas (px) in the 
chemical equilibrium case increases during expansion. This counterintuitive result 
appears due to the assumptions of entropy conservation and chemical equilibrium. 
The total number of particles in a fluid element decreases with proper time r due to 
mass effects in the hadron phase. Mass of particles contributes to the entropy den- 
sity significantly in low temperature and the proportionality between the entropy 
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Fig. 1. Mean transverse momentum for pions as a function of thermal freezeout 
temperature at 6 = 5 fm in y^SNN = 200 GeV Au+Au collisions. Solid (dashed) 
line shows a result with an assumption of an ideal, chemical equilibrium (chemically 
frozen) hadronic fluid. 

density and the number density violates in the temperature range under considera- 
tion. Then the total energy that a fluid element possesses initially is partly used for 
pdV work due to expansion as well as the case for chemically frozen fluids. As a fluid 
element expands, the remaining energy is distributed among the smaller number of 
particles in chemical equilibrium. These are the reasons why the different behav- 
ior of (pt) appears according to the assumption of chemical equilibrium/freezeout. 
Note that increase of (px) results in a "reheating" behavior of a fluid element: Tem- 
perature in the chemical equilibrium hadronic fluid drops more slowly than that in 
the chemically frozen fluid does as if there exists a sort of dissipation [^^- Under 
the chemical equilibrium assumption in ideal hydrodynamic simulations, increas- 
ing (pt) is commonly utilized so far to fix T^^ by fitting px slope. However, this 
is attained only by neglecting data of particle ratio. If particle ratios are fixed 
properly in hydrodynamic simulations to reproduce the data, pr slopes, especially 
for protons, are hardly reproduced [E] due to a lack of radial flow. The same is 
true for differential elliptic flow: dv2{pT) / dpr ~ V2/ (pt) is reproduced by canceling 
increasing behaviors of both V2 and (px) under chemical equilibrium assumption 
[ llOj. Agreement of the results from the CE model with pr spectra and V2{pt) 
data is reached only when one neglects particle ratio and the system keeps chem- 
ical equilibrium until kinetic freezeout. So, among three classes mentioned above. 
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Fig. 2. Pseudorapidity dependence of elliptic flow from hydro and hydro+hadronic 
cascade models are compared with data Dashed and dotted lines are results 
from ideal hydrodynamics with T*'' = 100 and 169 MeV respectively. Solid line is 
a result from a hydro+cascade hybrid model. 

the HC model turns out to be the only model which is able to reproduce particle 
ratio, pt spectra and V2{pt) data at once in a proper way. Therefore a picture of 
the dissipative hadronic corona together with the perfect fluid sQGP core, which is 
properly implemented in the HC model, is consistent with these experimental data 
observed at RHIC. 

3. 3D HYDRO AND HADRONIC CASCADE MODEL 

Establishment of a dynamical model for a detailed description of the whole space- 
time evolution is one of the main goals in the physics of relativistic heavy ion 
collisions. Dynamical models for description of the whole space-time evolution will 
be needed to draw the transport properties of the QGP from the experimental 
data. A dynamical framework has been developed [E] to describe three important 
aspects of relativistic heavy ion collisions, namely colour glass condensate (CGC) 
for collisions of two nuclei, ideal 3D hydrodynamics for space-time evolution of 
locally thermalised matter, and jet quenching for high pT non-thermalised matter. 
According to the discussion in the previous section, we incorporate a hadronic 
cascade model JAM [ 13, into our previous framework, the "CGC-|-hydro-|-jet" 
model [ ^1 , to establish a more realistic description of space-time evolution as 
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Fig. 3. Illustration of the approximately monotonia increase of absolute value of the 
shear viscosity with temperature. SYM, HRG, and wQGP represent, respectively, 
supersymmetric Yang-Mills model, hadronic resonance gas, and weakly coupled 
QGR 

discussed in the previous section. 

Figure 2 shows pseudorapidity dependences of V2 from this hybrid model and 
ideal 3D hydrodynamics with T**^ — 100 and 169 MeV. Here critical temperature 
and chemical freezeout temperature are taken as being Tc = T'^^ = 170 MeV in 
the hydrodynamic model. In the hybrid model, in which hydrodynamic description 
of the sQGP ideal fluid is followed by the kinetic theory of resonance gases, the 
switching temperature from a hydrodynamic description to a kinetic one is taken 
as Ts„ = 169 MeV. Ideal hydrodynamics with T'^ — 100 MeV which is so chosen 
to generate enough radial flow gives a trapezoidal shape of V2{ri) A large 

deviation between data [ 5 and the ideal hydrodynamic result is seen especially in 
forward/backward rapidity regions. However, just after hadronization (T**^ = 169 
MeV), an ideal sQGP fluid gives a triangle shape as shown by the dotted line in 
Fig. 12 (Note that the contribution from resonance decays dilutes slightly the elliptic 
flow generated in the QGP phase [El but that the shape is not so changed.) This 
means the large V2 in forward/backward rapidity regions in ideal hydrodynamic 
simulations with T*^ — 100 MeV is mainly generated in the perfect fluid of hadron 
phase. When hadronic rescattering effects are taken through the hadronic cascade 
model instead of perfect fluid description of the hadron phase, V2 is not so generated 
in the forward region due to the dissipation and, eventually, is consistent with the 
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Fig. 4. Illustration of the rapid variation of the dimensionless ratio of the shear 
viscosity, ^/(T), to the entropy density, s(T). Solid and Short dashed lines above 
Tc asymptotically merge wQGP result at high temperature region. We introduce a 
parameter which controls how fast the sQGP result reaches the wQGP one. Solid 
and short lines below Tc corresponds to sound velocity in the resonance gas model 
= 1/3 and 1/6 respectively. Long dashed line is rj/s = I/Att from SYM [ Igj. 

data. So the perfect fluid sQGP core and the dissipative hadronic corona picture 
works well also in the forward region [ 16 . Note that the eccentricity in the CGC 
initial condition is found to be larger than that in conventional participant scaling 
profiles. This is crucial to reproduce the V2 data in central collisions (3-15%) in 
our approach. For results from conventional parametrization based on Glauber 
approaches for the initial conditions, the deviation between ideal hydrodynamics [ 
121 and the data E El E] which was seen previously in peripheral collisions is 
interpreted only by the late hadronic viscosity [ .17) . On the other hand, the hybrid 
approach with CGC initial conditions overpredicts the f 2 data unless the early QGP 
viscosity is considered in peripheral collisions. This suggests that one needs a much 
better understanding of initial conditions [ 18 to draw the information about the 
transport properties of the sQGP from the elliptic flow data. 

4. Summary: What have we learned? 

At least in central collisions (up to ~ 20%), we can establish a new picture of 
space-time evolution of produced matter, namely "perfect fluid sQGP core and 
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dissipative hadronic corona" , from a careful comparison of hydrodynamic results 
with experimental data observed at RHIC. Note that, in semi-central to peripheral 
collisions, there is a room for existence of additional dissipative effects in the QGP 
phase from a recent quantitative analysis [E| of W2 data. This comes from presently 
unknown details of the initial stage in relativistic heavy ion collisions. Further 
theoretical efforts for understanding of the initial stage [ ^] are highly needed to 
obtain a more conclusive picture. In any case, a possibility of the picture, "perfect 
fluid sQGP core and dissipative hadronic corona" , is intriguing. 

What is the physics behind this picture? 77/s is known to be a good dimension- 
less measure (in natural unit h — ks ^ c ^ 1) to see the effect of viscosity, where 
r] is the shear viscosity and s is the entropy density. Figures 3 and 4 show possible 
scenarios for temperature dependence of rj and rj/s deduced from the discussion 
in the previous sections, rj/s must be small in the QGP phase, which might be 
comparable with a value l/An conjectured as the minimum one among any kinds of 
fluids [Ell and the perfect fluid assumption can be vahd. Whfle r]/s becomes huge 
in the hadron phase and the dissipation cannot be neglected. Shear viscosities of 
both phases are found to give 77 ^ 0.1 GeV/fm^ around Tc So shear viscosity 

itself is expected to increase with temperature monotonically. The "perfect fluid" 
property of the sQGP is thus not due to a sudden reduction of the viscosity at the 
critical temperature Tc, but to a sudden increase of the entropy density of QCD 
matter and is therefore an important signature of deconfinement. 
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